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Abstract — The most of the existing models assume a constant 
flow rate of river cross section during quality prediction. The 
application of such models in wide rivers can’t return the 
accurate results, which may not reliable for devising suitable 
mitigative measures. The river Brahmani is considerably much 
wide and many tributaries carry the effluents of different 
industries such as coal mines, thermal power plants and 
aluminium smelters etc with various pollutants, which have 
different characteristics in nature to degrade the river 
ecosystem. The present study has been conducted to simulate 
water quality along the stretches of Brahmani river in 
Angul-Talcher region of Orissa during winter season. The 
velocity variation over the cross section, changes in river 
dimensional, the impact of river course, and the percentage 
and the settling velocity of particulate maters duly 
incorporated in the present model. It has been observed that 
the model is effective and economic to study the water quality 
dispersion along the river stretches with 1.5 to 12.5% of error 
compare to the analyzed data. The present study may be useful 
for the concern industrialists and environmentalists to devise 
the suitable pollution control strategies. 

Index Terms — computation model, conservation 
parameters, flow rate, quality prediction, river water, velocity 
distributive function, 

I. INTRODUCTION 

The river water is being utilized for various purposes such 
as drinking, bathing, agriculture and industrial use. Most of 
the Industries are responsible for polluting the environment, 
particularly the nearby river systems by continuous discharge 
of the industrial effluents. Such Industries are always located 
on or nearby the bank of rivers. Therefore, it is very much 
warranted that a study on the dispersion characteristics of 
pollutant species along multi wastewater- outfall river system 
has to be considered in order to develop an indigenous and 
effective algorithm for water quality prediction, which would 
lead for designing suitable predictive system for better river 
quality management using appropriate mathematical models 
for simulating the water quality along the river stretches. 
Most of the existing software has many limitations. These 
models have been developed with many assumptions with the 
aim to simplify the governing equations into the standard 
form of differential equations, which can easily be solved 
through analytical methods. These models assume the 
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velocity over the cross section of the river system as constant. 
The assumption of constant velocity over the cross section 
would cause significant error in the measure of flow rate and 
the predicted water quality data would possess a less accuracy 
[Grantham, 1970; Thayer and Krutchkoff, 1967; 
Sundararajan et al. 2002a: Mohan and Sundararajan 2010b]. 
Therefore, a study has been carried out for predicting 
non-conservative parameters along the 

multi-industrial-wastewater-outfall river stretch using the 
basic Streeter and Phelps model [Streeter and Phelps, 1925], 
which was improved by incorporating the ratio of setteable 
bio-flocculated particulate matters [Bhargava, 1989] and 
later it was further improved duly incorporating velocity 
variation over the cross section of the river system and river 
dimensional variations adopting various steps, techniques 
and tools of system analysis [Sundararajan et al. 2002a]. 

In this research paper, an appropriate method of predicting 
surface water quality using various steps, techniques and 
tools of system analysis has been presented along with a case 
study carried out for predicting water quality along the 
stretches of Brahmani River in Angul-Talcher Region of 
Orissa. 

II. MATERIALS AND METHODS 

A. Study Area 

The study area is bounded by latitudes 20° 37’ 00” N to 
21° 10’ 00” N and longitudes 84° 53’ 00” E to 85° 28’ 00” E 
and situated at an average height of 139 m above mean sea 
level. The river catchment is characterized by Precambrian 
granites, gneisses and schists of Eastern Ghats with local 
intrusive and volcanic lithologies; lime stone, sand stone and 
shales of the Gondwanas (Konhauser et al. 1997; Panda et al. 
2006; Sundaray et al. 2006; Sundaray 2009). River 
Brahmani enters the Angul district through Rengali reservoir 
and passes through Angul-Talcher-Meramundali region 
(from northwest to southeast) of Orissa. Finally Brahmani 
river drains into Bay of Bengal after flowing through Jajpur 
& Kendrapara district of Orissa. The drainage network is 
controlled by river Brahmani. It is a perennial river and has 
innumerable large and small tributaries. It meets with its four 
tributaries such as Tikira, Singhara and Nandira in 
Angul-Talcher region of Angul district while Kisinda meets 
the Brahmani river at Meramundali in Denkanal district. 
Water of Brahmani and its tributaries cater to the 
industrial/domestic need of this fast growing industrial 
complex and in return they release thousands of gallons of 
wastewater to the river and tributaries. The drainage pattern 
is controlled by the River Brahmani along with its tributaries. 
River Brahmani plays as a major source of water, as well as 
sink for discharging effluents from the existing industries. 
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Fig.l. Geographical representation of study area 

river stretches. Further, the Bramani river has many 
tributaries, which carry the pollutants and toxic elements 


B. Sampling and analysis 
In order to have a fairly comprehensive picture of the 
intensity of pollution of the river systems, winter season has 
been considered. For the river dimension can be measured 
more accurately in winter season than other seasons. In the 
monsoon season, the pseudo data are inputted as river 
dimensions due to excess water flow while the flow becomes 
very lean due to the high evaporation rate and almost rainless 
days in the summer season. The water samples were collected 
from all five tributaries, just before meeting into Brahmani 
river. The Samples were taken from 10 to 15cm below the 
water surface using previously cleaned polyethylene 
container as per standard procedures (APHA, 2005). Each 
analysis was carried out in triplicate and the mean value was 
taken for consideration. The physical parameters of river 
such as width, depth and stream velocity the tributaries as 
well as the main river were considered for the estimation of 
accurate and rational flow rate. The entire river course 
between the reference point at Anantapur and the end point 
at Kisinda Jhor has been divided into six stretches by means 
of the mixing points of five tributaries i.e., the Tikira, the 
Singhara, the Bangaru, the Nandira and the Kisinda, and 
further the model has been applied to predict the water 
quality along these river stretches. 


from different sources such as coal mines, TPPs, iron & steel, 
aluminium industries and townships etc, The selected model 
incorporates the velocity variation over the cross section of 
the river, dimensions of the river stretches such as width, 
depth and stream velocities. It also takes into account the 
physiochemical characteristics of tributaries that blend with 
the main river. 

The model has been developed based on fundamental truth 
that the concentration of the non-conservative water quality 
parameter at time ‘t’ is directly proportional to the rate of 
change in the initial concentration. The measure of BOD 
contributed by bio-flocculated particulate matters has been 
accounted in this study. The model has been derived taking 
the settling velocity of the settleable bio-flocculated 
particulate matter. The flow rate is one of the most important 
input parameters in river quality prediction. No model is 
appropriate to assess the accurate flow rate of the river 
system. The equations/software are used in assess the flow 
rate with very less percent of error from data obtained in the 
field survey. Therefore, an appropriate model has been 
chosen for computing accurate and rational flow rates in 
order to avoid error in prediction. 


C. Prediction modeling 

The climate condition of the study area does not vary 
dramatically during winter. Therefore the effect of solar 
radiation is assumed to be negligible in this area. The quality 
of water is significantly controlled by the river water flow as 
well as varying nature of river dimensions such as river 
width, depth and river course, which considerably alter the 
stream velocity. As a result, the dispersion scenario of water 
pollutants in the river stretches considerably vary along the 


1) Velocity distributive function 
The velocity distributive function in the cross section 
of the water flow region in laminar flow condition has been 
derived as 


V(x, y,z) = V 
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In case of the river whose cross section is 
non-symmetric subject to depth axis, the velocity of the 
stream flow may be assumed by the following function: 


V(x, y,z) = 
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1 ) Computation of flow rate 

If u (x,y) is the average stream velocity of the finite 
element whose center is (x,y), then the flow rate of river or 
stream can be computed as follows: 


D 

Q = 2 Z 

y = Sy/2 
step = Sy 


Z [S(x,y)(5x.<%0] 
x = 0 
step = <Sx 

(3) 


Estimation of dispersion coefficient 

The extent to which a stream can be approximated by a one or 
two-dimensional model also influences how well dispersion 
coefficients can be estimated. Furthermore, the degree of 
averaging in space and time in solving the 
affective-dispersive equation can influence the magnitude of 
the dispersion coefficient shown by calibration. In some 
cases, numerical dispersion can be comparable, to actual 
dispersion. When this occurs, the dispersion is specified as 
zero (i.e., Dx=0), or better still as the difference between the 
actual and the numerical dispersion. Dispersion coefficients 
can be estimated for every river and such estimation may 
quantify the ability of the river to assimilate waste. Later 
dispersion studies were conducted to show that dispersion 
was not an important part of the steady-state waste assimilate 
capacity. More recently, interest seems to have returned to 
the dispersive capability of a stream as concern have 
heightened regarding the mixing of hazardous waste spills 
and as efforts have intensified in the modelling of dynamic 
water quality conditions. The dispersion coefficient (k) for 
the river can be estimated using the following mathematical 
formula. 


k 





( 6 ) 


In case of non-symmetric cross sectional river system 
with the extremities of top layer width (a G , 0) and (b 0 , 0), the 
flowrate can be computed as follows: 

The flow rates of the river in the positive and 
negative direction of width axis have respectively been 
computed using the following Equations (4) and (5). 

Thus the total flow rate of the stream over the 
non-symmetric cross sectional river system can be obtained 
as Q=Q + +Q . 

When the finite element thickness 5x and 5y are tending to 
zero, Q tends to accurate. 

f \ 


where u is average stream velocity, tl and t2 are the mean 
times of passage for the dye cloud to move past stations 1 and 
2, and a2tl and a2t2 are the variance of the time versus 
concentration curves. 

Influence of settleable BOD component 
The fraction of settleable BOD component contributes while 
it is transported along the river stretch up to some distance 
depending upon the settling velocity of the settleable 
components in dynamic condition. If the settling velocity is 
significantly less, then the distance along the downstream 
stretch will be increasing. Bhargava has attempted to 
incorporate the contribution of settleable BOD component in 
order to predict the concentration of part of BOD 
concentration, Cs with respect to the fraction of settleable 
BOD component, p. 
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where vs is the settling velocity of the component, Co is the 
initial concentration, D is the average, depth of the river and 
t is the time travelled by the components along the river 
stretch in the downstream. Further, Bhargava [1989] has also 
derived the other component of BOD prediction, resulted by 
the non-settleable BOD component as follows: 


C =(l-p)C e" kt 
ns o 


( 8 ) 
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The above Equations (7) and (8) are added up to get the 
combined effect of settleable as well as non-settleable BOD 
components in BOD prediction along the river stretch. In the 
above attempt, the influence of stream velocity on the settling 
velocity of the setteable particulate matter in dynamic river 
system has not been accounted. The model has been further 
improved by modifying as a 3D model, taking the variation of 
velocity and depth over the cross section of river system as 
follows. 


Three-dimensional water quality prediction model 
Bhargava (1989) have been derived to predict the water 
quality along the river stretch assuming the river stretch as 
one-dimensional and the depth and width of the river are 
constant. This model can be easily modified as a 
three-dimensional if the cross section of the river is divided 
into finite number of volumetric element in order. By using 
the velocity distributive function, the velocity of each 
volumetric element can be predicted. Since the length of each 
stretch is known, the varying time length taken by each 
volumetric element to travel the constant length of the stretch 
can be computed. Here it can be noted that the number of 
elements over the cross section will be constant for each 
stretch of the river. It is assumed that the dimensions of the 
volumetric element corresponding to either width or depth 


C(x, y, z) = p\ 


c o v(x, y,z ) 
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-k v(x,y,z) 






axis increases when there is an increase in the respective 
dimensions of width or depth. In this case, it is necessary to 
assess the flow rate of water passing the cross section of each 
volumetric element. Since such tasks involve number of 
calculations and incorporation of depth and width variation 
of stretches, the river stretch can be sufficiently divided into 
finite number of segments. However, the computer program 
has to be designed to incorporate all the data related to the 
river dimensions as well as the in-taking and discharging 
rate of water and effluent and the physico-chemical 
characteristics of the effluent of the point source mixing with 
the main river system. If there is no such source, except the 
river dimensions all other parameters can be taken as zero. 
The 3D-water quality prediction model for predicting the 
non-conservative parameters, whose initial concentration 
and maximum stream velocity just after the critical mixing 
distance are co and vo respectively at the location x=0, y=0, 
z=0 has been derived adding the Equations (7) and (8) in 
order to have the combined effect of settleable and 
non-setteable particulate matters. Thus, the 3D water quality 
prediction model may be obtained duely incorporating the 
effect of settling velocity of settleable bio-floculated 
particulate matters as follows: 


v(x,y,z) 


kz 


+ ( 1-/0 


c 0 v(x, y,z ) 


v(x, y,z) 


where C(x,y,z) is the concentration of the non-conservative 
parameter in the water transported through the volumetric 
element whose centre is the point P(x,y,z) and the velocity is 
v(x,y,z). p is the fraction of the bio-flocculated particulate 
matters, v’s is its settling velocity in static condition of the 
river, D is the maximum depth of the river over the width but 
average over the length of the stretch, k is the dispersion 
coefficient of the pollutant species and k’ is the dispersability 
coefficient of the settleable bio-flocculated particulate 
matters. First, the computer program input the data related to 
the reference point of the river. Successively, the subroutine 


starts for n times where n is the number of outfalls considered 
for the prediction studies and the program inputs the 
dimensions of the river and the physico-chemical parameters 
related to respective outfalls for the respective number of 
subroutines in order to predict the concentration of the water 
quality along the successive stretches. The flowchart of the 
computer program for water quality prediction along the 
river stretch joining the points, where the stream velocity is 
more than the velocities over the cross section of the river 
system has been presented in Figure 2. 



Fig. 2 Schematic diagram for river stretch along with their tributaries 
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R = River flow rate of at reference point 
C = Water quality of river at reference point 

X = Distance from last outfall (Xl=2 km, X2=5km, X3=19.5km, X4=36.5km, X5= 49km) 

S =Name of Outfalls/Streams (S0= Reference point, Sl=Tikira Stream, S2= Singhara Stream, S3=Bangaru Stream, 
S4=Nandira Stream, S5= Kisinda Stream 
W = Water quantity drawn from outfalls 
q = Flow rate of each outfall and 
b = Water quality of outfall. 


Table. ! Details of input data along with tributaries of Brahmani river 


Code 

^Parameters 

RiVcts 

Longitude 

Latitude 

Distance 
from last 
out fall 
(km) 

Total Width 
(W) 

Maximum 
Depth (Dm) 

Maximum 

Velocity 

(Vm) 

ao 

bo 

Flow Rate 
(m3/sec) 

SO 

Reference Point 
(near Ananta pur) 

85007’19” 

21006W’ 

0.0 

55.0 m 

1.06 m 

0.14 m/s 

-30 

25 

2.57 

SI 

Tikira river 

85005’07” 

21005’40” 

2.0 

35.0 m 

1.02 m 

0.25 m/s 

-15 

20 

1.15 

S2 

Singhara river 

85006’50” 

21004W’ 

3.0 

14.0 m 

0.73 m 

0.13 m/s 

-4 

10 

0.352 

S3 

Bangaru river 

85014W’ 

21000’40” 

14.5 

10.0 m 

0.61 m 

0.17 m/s 

-0.5 

9.5 

0.274 

S4 

Nandira river 

85015’35” 

20053W’ 

17.0 

20.0 m 

0.76 m 

0.4 m/s 

-10 

10 

1.613 

S5 

Kisinda river 

85020’30” 

20048’30” 

12.5 

8.5 m 

0.67 m 

0.14 m/s 

-5 

3.5 

0.174 


Table 2. Water Quality of various Tributaries of Brahmani river in Winter Season 


SI. No 

Parameters 

Ref. Point 

Tikira 

Singhara 

Bangaru 

Nandira 

Kisinda 

(SO) 

(SI) 

(S2) 

(S3) 

(S4) 

(S5) 


TDS (mg/1) 

174 

206 

411 

594 

288 

341 


ss (mg/1) 

41 

74 

156 

216 

56 

97 


Fluoride (mg/1) 

0.1 

0.1 

0.9 

1.2 

1.3 

1.8 


Chloride (mg/1) 

7.1 

14.2 

67.2 

46.2 

10.7 

28.4 


Sulphate (mg/1) 

40 

28 

52.8 

62.9 

25.0 

42.5 


DO (mg/1) 

7.0 

7.8 

4.9 

4.3 

6.7 

4.2 


BOD (mg/1) 

1.8 

1.9 

2.7 

4.9 

1.3 

1.7 


COD (mg/1) 

27.0 

44.6 

107.3 

126.5 

66.8 

86.3 


III. RESULTS AND DISCUSSION 

As the concentration of non-conservative parameters change 
due to dispersion phenomena, the concentration of such 
parameters varies with respect to the ‘time’ or ‘distance’. 
Using the mathematical formula, which has been discussed 
elsewhere in the thesis, the dispersion coefficients have been 
estimated and used for the prediction of non-conservative 
parameters as follows: k = -2.0 x 10-5 for BOD, COD, TSS, 
S04, Fluoride and Chloride, k = -2.0 x 10-7 and k = 1.4 x 
10-7 DO and TDS. The calibrated model was used to predict 
the water quality with an independent set of data as a part of 
validation exercise. Results of model predictions were fairly 
good, and performance of the model was further confirmed 
through statistical evaluation of the results. 


A. Biological Oxygen Demand 

A reference point two km away from the confluence point of 
Tikira tributaries with Brahmani river was considered at 
Anantapur. Initially the concentration of BOD was 1.8 mg/1 
and it was decreasing in the downstream into 1.1 mg/1 and 
then mixed with the Tikira tributary. The concentration after 
the critical mixing distance increases to 1.3mg/l, which has 
been predicted to be 1.43 mg/1 with acceptable error of 10% 
using the present model. Again the concentration of BOD 
decreases gradually into 0.78 mg/1 after travelling of 3 km 
distance and then mixed with Singhara tributary whose BOD 
concentration was 2.7 mg/1. The concentration of BOD in 
Brahmani hikes into 1.45 mg/1 after the critical mixing 
distance, whereas the predicted data shows to be 1.41mg/l 
with error 2.8% which is acceptable. Now the BOD 
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concentration decreases significantly into 0.16 mg/1 after 
travelling of 14.5 km. The concentration of BOD in 
Brahmani river goes up into 0.53 mg/1 after the mixing 
distance when Bangaru tributary blended with the main 
river. In the case the predicted concentration is 0.55 mg/1 
with error of 3.6%. The concentration of BOD in Brahmani 
river after following of 17.5 km significantly decreases and 
blends up with Nandira tributary whose BOD concentration 
is 2.4 mg/1. After well-mixing, the concentration of BOD in 
Brahmani becomes 0.75 mg/1, which has been predicted as 
0.7 mg/1 with error of 6.6%. The concentration gradually 
decreases and mingles with Brahmani and results into the 
concentration of 0.23 mg/1, which has been predicted as 0.22 
with error of 3.4%. The concentration of BOD was predicted 
for five kilometers after the confluence point of Kisinda and 
Brahmani and it has been observed from the result that the 
quality of water has been significantly improved. In all the 
above cases, the error percentage varies from 2.8 to 10, which 
may be acceptable for practical purpose of river quality 
management. 



SO SI S2 S3 S4 S5 0 

Distance (m) 

Fig. 3 Prediction of BOD concentration within river stretch 
B. Dissolved Oxygen 

Initially the concentration of DO at the reference point of 
river was 6.9 mg/1. In actual analysis of DO, the 
concentrations of DO in the downstream of Brahmani river 
are: 7.3 mg/1 just after the confluence of Tikira and 
Brahmani, 7.2 mg/1 just after the confluence of Singhara, 7. 1 
mg/1 just after the confluence of Bangaru, 6.9 mg/1 just after 
the confluence of Nandira and 7.1 mg/1 just after the 
confluence of Kisinda, whereas the predicted concentration 
of DO are 7.3 mg/1, 7.1 mg/1, 7 mg/1, 7 mg/1 and 6.9 mg/1 at 
the respective confluence points of the tributary with the 
main river. In all the cases, the error varies from 0 to 2.8%, 
which may be acceptable. 



Fig. 4 Prediction of DO concentration within river stretch 


Total Dissolved Solids 

The actual concentration of TDS, near the reference point, 
was 194 mg/1 and it was decreasing in the downstream and 
then it mixed with the Tikira tributary. The concentration of 
river becomes 199 mg/1. After the mixing of Singhara and 
Bangaru tributaries the predicted TDS concentration of river 
goes up to 216 and 210 mg/1 while the actual data was 209 
and 214 mg/1 respectively. After travelling of 17.5 km 
distance the concentration was reduced to 205 mg/1 and then 
mixed with Nandira (233 mg/1), the predicted TDS 
concentration of river becomes 211 mg/1 while the actual 
concentration of the river water was 202 mg/1 with 1.8% 
error. The concentration has been predicted as 208 mg/1 after 
the mixing of Kisinda whereas the actual concentration was 
211 mg/1 with error of 4.4%. Predicted model results indicate 
that, simulated TDS results were relatively in best agreement 



Fig. 5 Prediction of TDS concentration within river stretch 


C. Total Suspended Solids 

TSS value near the reference point was 41 mg/1 and 
decreased into 31.8 mg/ in downstream. After the mixing of 
Tikira the maximum predicted concentration of TSS in river 
water was 49 mg/1 whereas the actual was 45.5 mg/1 with 
error of 7%. Similarly the TSS concentration was predicted 
as 43.7 mg/1 after the mixing of Singhara while the actual 
data was 46.5 mg/1 with 6% error. The TSS concentration 
were decreasing (6.4 mg/1) with the distance of 14.5 Km, 
then mixed with Bangaru and the concentration was raised 
up to 1 1.2 mg/1 with maximum error of 9%. Again after the 
travelling of long distance (17 km), then the river water with 
the TSS concentration of 1.2 mg/1 mixed to the Nandira and 
concentration of TSS was going to increased up to 16.2 mg/1 
where as the actual concentration was 14 mg/1 with error of 
7%. After the mixing of Tikira the TSS concentration was 
predicted as 6.5 mg/1 while actual concentration was 7. 1 mg/1 
with 7.6 % error. 



Fig. 6 Prediction of TSS concentration within river stretch 
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D. Chemical Oxygen Demand 

In the case of COD, initially the concentration of COD was 
57 mg/1 and it was decreasing in the downstream into 44.2 
mg/1 and then mixed with the Tikira tributary. The 
concentration after the critical mixing distance increases to 
51mg/l, which has been predicted to be 55 mg/1 with 
acceptable error of 7.8%. Again the concentration of COD 
decreases gradually into 37.3 mg/1 after travelling of 3 km 
distance and then mixed with Singhara tributary whose COD 
concentration was 127 mg/1. The concentration of COD in 
Brahmani hikes into 43 mg/1 after the critical mixing 
distance, whereas the predicted data shows to be 44.5 mg/1 
with error 3.4% which is acceptable. Now the COD 
concentration decreases significantly into 6.3 mg/1 after 
travelling of 14.5 km. The concentration of COD in 
Brahmani river goes up to 9 mg/1 after the mixing distance 
when Bangaru tributary blended with the main river, while 
the predicted concentration was 8.4 mg/1 with error of 6.6%. 
The concentration of COD in Brahmani river was 
significantly decreased after following of 17.5 km and blends 
up with Nandira tributary whose COD concentration is 76.2 
mg/1. After well-mixing, the concentration of COD in 
Brahmani becomes 19.5 mg/1, which has been predicted as 
18 mg/1 with error of 7.7%. The COD concentration 
gradually decreases and mingles with Kisinda and results 
into the concentration of 6.2 mg/1, which has been predicted 
as 5.6 with error of 9.6%. The concentration of COD was 
predicted for five kilometres after the confluence point of 
Kisinda and Brahmani and it has been observed from the 
result that the quality of water has been significantly 
improved. In all the above cases, the error percentage varies 
from 3.4 to 9.6 which maybe acceptable for practical purpose 
of river quality management. 


Predicted —■—Analysed 



Distance (in) 

Fig. 7 Prediction of COD concentration within river stretch 
E. Chloride 

Chloride concentration at the reference point of river was 
14.2 mg/1. In actual analysis of chloride, the actual 
concentrations of chloride in the downstream of Brahmani 
river are: 15.5 mg/1 (after the confluence of Tikira and 
Brahmani), 15.2 mg/1 (after the confluence of Singhara and 
Brahmani), 3.6 mg/1 (after the confluence of Bangaru and 
Brahmani), 2.8 mg/1 (after the confluence of Nandira and 
Brahmani) and 1.35 (after the confluence of Kisinda and 
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Brahmani), whereas the predicted concentration of chloride 
are 16.3mg/l, 14.3mg/l, 3.3 mg/1, 2.65 mg/1 and 1.3 mg/1 
respectively. In all the cases, the error varies from 6 to 9% 
which may be acceptable. 



Distance (m) 

Fig. 8 Prediction of Chloride concentration within river 
stretch 

F. Sulphate 

The actual concentration of sulphate near the reference point 
was 28.3 mg/1 and it was decreasing in the downstream and 
then it mixed with the Tikira tributary. The sulphate 
concentration of river was 30.9 mg/1 while the predicted 
concentration becomes 32.2 mg/1 with error of 4%. The 
concentration of sulphate decrease gradually in to 21.8 mg/1 
and after the mixing of Singhara and Bangaru tributaries the 
predicted sulphate concentration of river goes up to 26 and 7 
mg/1 while the actual concentration was 28 and 8.5 mg/1 
respectively. After travelling of 17.5 km distance the 
concentration was reduced to 0.75 mg/1 and then mixed with 
Nandira (45 mg/1). The actual sulphate concentration was 
13.5 mg/1 where as the predicted concentration of river 
becomes 12.8 mg/1 1 with 5.5% error. The concentration has 
been predicted as 4.4 mg/1 after the mixing of Kisinda 
whereas the actual concentration was 4.8 mg/1 with error of 
10 %. 



Fig. 9 Prediction of Sulphate concentration within river 
stretch 
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G. Fluoride 

Fluoride’s concentration near the reference point was 0.9 
mg/1 and decreased into 0.7 mg/1 in downstream. After the 
mixing of Tikira the actual concentration of fluoride in river 
water was 0.74 mg/1 whereas the predicted was 0.72 mg/1 
with error of 4%. Similarly the fluoride concentration was 
predicted as 0.55mg/l after the mixing of Singhara while the 
actual data was 0.5mg/l with 9.3% error. The fluoride 
concentration were decreasing (0.08 mg/1) with the distance 
of 14.5 km, then mixed with Bangaru and the actual 
concentration was raised up to 0.24 mg/1 while the predicted 
concentration was 0.22 with error of 9%. Again after the 
travelling the distance of 17 km, the river water with fluoride 
concentration of 1.2 mg/1 mixed to the Nandira then the 
concentration was increased up to 0.37mg/l where as the 
actual concentration was 0.34mg/l with error of 8%. After 
the mixing of Tikira the fluoride concentration was 0. 13mg/l 
while prediction concentration was 0.14 mg/1 with 7.6 % 
error. 



Distance (in) 

Fig. 10 Prediction of Fluoride concentration within river 
stretch 

IV. CONCLUSION 

The model has been applied for predicting river water quality 
during winter only as turbulence flow takes place during 
monsoon and lean flow during summer season. The selected 
model is only applicable while laminar flow takes place. 
Water quality prediction during lean period in Bramani river 
is somewhat complex as diversified flows take place in the 
wide river bottom. 

The model is effective in predicting the non-conservative 
parameters with errors varying from 0 to 10%. In the cases of 
predicting the non-conservative parameters, the error 
percentage varies from 0 to 10%, which is acceptable in 
predicting water quality in mining and its allied industrial 
complexes [Sundararajan, 2004]. Here the error percentage 
may again decrease if the non-point sources are considered in 
the present study as the present study has been carried out 
only considering the point sources. 
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